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ABSTRACT 

Transonic discs with accretion rates relevant to intrinsically bright Galactic X-ray 
sources (L rj 10 38 -10 39 ergs -1 ) exhibit a time dependent cyclic behaviour due to the 
onset of a thermal instability driven by radiation pressure. In this paper we calculate 
radiation spectra emitted from thermally-unstable discs to provide detailed theoretical 
predictions for observationally relevant quantities. The emergent spectrum has been 
obtained by solving self-consistently the vertical structure and radiative transfer in 
the disc atmosphere. We focus on four particular stages of the disc evolution, the 
maximal evacuation stage and three intermediate stages during the replenishment 
phase. The disc is found to undergo rather dramatic spectral changes during the 
evolution, emitting mainly in the 1-10 keV band during outburst and in the 0.1-1 
keV band off-outburst. Local spectra, although different in shape from a blackbody 
at the disc effective temperature, may be characterized in terms of a hardening factor 
/. We have found that / is rather constant both in radius and in time, with a typical 
value ~ 1.65. 



Key words: 

transfer 



accretion, accretion discs - black hole physics - instabilities - radiative 



1 INTRODUCTION 

Observations of Galactic black hole candidates (BHCs) and 
active galactic nuclei (AGNs) have gone a long way in ob- 
taining very detailed information on these sources. The anal- 
ysis of X-ray data revealed the existence of at least four 
distinct spectral states in BHCs, defined on the basis of 
the spe ctral shape and the flux level (see e.g. Tanaka & 
Lewin (1995); Narayan, Mahadevan & Quataert (1999) for 
reviews). In the high (soft) state the spectrum is dominated 
by a soft thermal component, commonly attributed to a 
standard accretion disc. In the low (hard) state a power-law 
high energy tail, probably originated by the comptonization 
of soft photons by a hotter medium, dominates the spec- 
trum. Usually, as the flux decreases, the power-law compo- 
nent becomes more and more pronounced, while the thermal 
one weakens and eventually disappears. Optical/UV/X-ray 
AGN spectra consist of three main components: a big blue 
bump, a soft X-ray excess, a hard power-law tail with a re- 
flectio n hum p and a high-energy cut-off (see e.g. Zdziarski 
et al. (1996) for a review). As it was realized long ago, by 



the energy shift) is corroborated by a nu mber of observa- 
tional evidences (see e .g. Gr ebenev et al. (1993); Grebenev, 
Sunyaev & Pavlinsky (1997) for the properties of power-law 



tails and Zdziarski, Lubihski & Smith (1999) for the spectral 



shifting a typical AGN spectrum to higher energies by a fac- 
tor ~ 100, such that the big blue bump falls in the 1-10 keV 
band, one obtains a X-ray spectrum which closely resembles 
that of a typical BHC in an intermediate state. The ac- 
tual similarity between these two classes of sources (modulo 



index-Compton reflection hump correlation). 

In the last few years, a great wealth of information on 
the physical conditions in BHCs has come from the study 
of their variability. Besides transitions among different spec- 
tral states (on time-scales w 1 hour-1 month), time analysis 
revealed much more rapid variations, ranging from very low 
(mHz) to high-frequency (100 Hz) QPOs. In some sources 
(e.g. GRS 1915+105) phase lags between photons in differ- 
ent X-ray bands have been detect ed, w ith typical va lues i n 
the range 10-1000 ms (Cui et al. ( |1999| ); Reig et al. (|2000| ); 
Cui, Zhang & Chen ( fcOOOp ; Pottschmidt et al. ( |200C| )). If 
the characteristic times scale linearly with the black hole 
mass, as one naively expects, the shortest variability in 
AGNs should occur on ~ hours-days. Up to now, no def- 
inite evidence for QPOs nor phase lags has been reported 
in AGNs, although multi-wavelength observations of some 
active galaxies seem to indicate a correlated variabil ity be - 
tween the opti cal- UV and X-ray flux es (C lavel et al. ( 1992 ); 
Edelson et al. jl996| ); Nandra et al. ( |l99g| )). Time delays of 
a few days between the optical-UV and X-ray bands could 
be detectable, providing a powerful tool to probe the geom- 
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etry of the accretion flow in t he nu clear regions, as recently 
suggested by Zampieri et al. (2000). 



The analogies between the spectral properties in the 
optical/UV/X-ray bands and, possibly, the variability pat- 
terns in BHCs and AGNs suggest that in both classes of 
sources the same basic mechanism might be responsible for 
the observed emission. Much effort was devoted in recent 
years to envisage a "unified" scenario that has been par- 
tially successful in explaining the different spectral states in 
BHCs. The current picture is based on a two-component ac- 
cretion flow: a standard disc and a hotter phase which could 
be advecti on-do minated (see e.g. Narayan, Mahadevan & 
Quataert (f999)). The different states are then produced by 



the varying importance of the two components. However, 
specific models do not provide a consistent description for 
the structure and evolution of the accretion flow. 

It has been known since a long time that the inner re- 
gions of a standard disc are thermall y and viscously unsta- 
ble (e.g. Pringle, Rees & Pacholc zyk ( 1973); Shapiro, Light- 
man & Eardely fll976| ); Ichimaru ( p.977[ )). This suggested the 
intriguing possibility that the observed spectral transitions 
may be triggered by such instabilities. A definite answer 
should be sought in terms of a time-dependent hydrodymi- 
cal appr oach, such as that adopted by Ho nma, M atsum oto 
& Kato (|1991[) and Szuszkiewicz & Miller (EgOT), fll998|). In 



particular, the latter authors have calculated the non-linear 
(global) evolution of transonic accretion discs produced by 
the onset of a thermal instability driven by radiation pres- 
sure and have shown that the evolution follows a cyclic pat- 
tern, passing through successive stages of evacuation and 
replenishm ent of the inner disc region (see also Takeuchi & 
Mineshige (|I998|)). 



Ultimately, in order to get detailed theoretical predic- 
tions for the observed quantities, it is necessary to com- 
pute the spectra emitted during the instability cycle. Model 
spectra should properly take into account for the effects 
of the disc atmosphere, through which energy deposition 
is likely to take place at different heights. Models of dis c 
atmospheres have bee n alrea dy presented (LaDous (1989); 
Shimura & Takahara (1993)) and some thought has been 
given to the prob lem o f energy deposition (Shaviv & Wehrse 
(1986); Hubeny ( 199C )) . The first thorough investigation of 



the vertical structure and emission spectrum of a thin accre- 
tion disc has been presented by Shimura & Takahara ( p.995 ), 
who found that, for intermediate accretion rates, local spec- 
tra, emitted from each disc annulus, can be approximated by 
a Planckian with an average hardening factor / ~ 1.7. Inves- 
tigations of the spectra emitted by steady-state, dissipative 
accretio n dis cs in AGNs have been carried out by Sin cell & 
Krolik (|l998j) (standard discs) and Wang et al. (|l999|) (slim 



discs). Recently, Hubeny and co-workers (( 2000[ ) and refer- 
ences therein) reported on the progress in their long term 
programme to construct detailed spectra from stationary, 
thin accretion discs in AGNs. 

In this paper we compute spectra under the slim-disc 
approximation, which is appropriate for modeling intrin- 
sically bright sources and allows us to follow the non- 
stationary evolution of a disc subject to a thermal instability. 
We will assume that energy deposition takes place entirely 
in the disc. More complex situations in which energy depo- 
sition occurs also in the atmosphere will be considered in a 
follow-up paper. This work is part of an ongoing research 



programme dedicated to the analysis of hydrodynamical in- 
stabilities in slim accretion discs and to the modeling of the 
spectral and variability properties of BHCs, in particular of 
the microquasar GRS 1915+105. In this paper we compute 
the spectra that emerge from slim discs during their insta- 
bility cycle and discuss their possible relevance for the soft 
spectral component of GRS 1915+105. The paper is orga- 
nized as follows. In Section ^| we describe the structure and 
evolution of thermally unstable transonic disc, concentrat- 
ing on those properties relevant to spectral calculations. In 
Section [] our approach to the solution of radiative transfer 
in the disc atmosphere is outlined. Computed spectra are 
presented in Section ^| Discussion and conclusions follow in 
Section |^. 



2 RADIAL STRUCTURE OF A THERMALLY 
UNSTABLE TRANSONIC ACCRETION 
DISC 

Thermally unstable discs may undergo a limit-cycle be- 
haviour with successive evacuation and refilling of their cen- 
tral parts. By means of a full time-dependent hydrodynam- 
ical simulation, Szuszkiewicz & Miller (1998) explored in 
detail a model with black hole mass M — IOMq , initial ac- 
cretion rate M = 0.96 (in Eddington units) and viscosity 



parameter a = 0.1. In this particular model the cycle has 
a duration of about 780 s and can be divided into two dis- 
tinct parts: a very short (~ 20 s) outburst phase, during 
which the inner part of the disc is evacuated, and a quies- 
cent stage (corresponding to the slow refilling) lasting for 
~ 760 s. We have chosen four particular snapshots during 
the disc evolution for which to calculate emergent spectra. 
They correspond to the end of the refilling phase, just before 
the outburst starts (model 1), the peak of the outburst (the 
maximal evacuation stage, model 2), the beginning (model 
3), and half way through the refilling phase (model 4). The 
theoretical bolometric light curve together with the posi- 
tions of the four chosen snapshots is shown in Fig. [lj In this 
section we discuss the radial structure of the disc for each 
of the four stages. 

The electron temperature distributions for models 1- 
4 are shown in Fig. 0; unless stated otherwise, hereafter 
all lengths are in units of the Schwarzschild radius, re = 
2GM/c 2 . Before the instability sets in (model 1) the temper- 
ature profile is that given by the solid line, labeled 1 in Fig. 
|j| At the onset of the instability, two density waves are sent 
out from the locally-unstable region, one moving inwards 
and the other outwards. At the same time the temperature 
rises significantly in the unstable region. As the outgoing 
wave progresses the temperature peak is reduced (model 2) 
but remains still above that of model 1. Another effect of the 
propagating wave is to push a significant amount of matter 
outwards leaving behind an underdense region. This can be 
seen in Fig. |^ where the volume density is plotted against ra- 
dius. The outgoing wave heats the material through which it 
passes, causing the perturbed part of the disc to swell up, as 
shown in Fig. ||. During the evacuation phase the disc tem- 
perature is high. However, once the wavefront has moved 
beyond the linearly unstable region, the disc cannot remain 
at such a high temperature for a very long time. Eventually 
the front weakens and the temperature falls dramatically 
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Figure 1. The bolomotric light curve with the positions of four 
particular evolutionary stages described in the text. 
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Figure 2. The radial profile of the midplane electron temperature 
for the four particular evolutionary stages discussed in the text. 



(model 3). After this, the underdense region begins to fill 
up on a viscous time-scale. This increase in density is asso- 
ciated with a progressive rise in the temperature. 

A quantity of immediate interest for spectral calcula- 
tions is the effective optical thickness of the disc, r e // = 
b~a(T"es + Ta)} 1 / 2 ; here r ea and r a are the scattering and ab- 
sorption depths respectively. As shown in Fig [|, during the 
the maximal evacuation stage (model 2) the effective opti- 
cal depth drops slightly below unity in the disc inner region. 
However, r e ff is sufficiently large to safely assume that the 
radiation field in the disc is in LTE. Then, close to the equa- 
torial plane the radiation intensity approaches a blackbody. 
In traversing the disc atmospheric layers radiation goes out 
of equilibrium and scattering effects are likely to change the 
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Figure 3. Same as in Fig. ^ for the volume density. 




Figure 4. Same as in Fig. |2| for the disc half-thickness 



spectrum, as we will discuss in more detail in Section ]^ 
The total flux distribution for each stage of the evolution 
is shown in Fig. ^. Neglecting relativistic effects, the total 
luminosity of the disc in the four different stages Ldisc can 
be computed integrating the flux over the disc surface (see 
Table |l|). It should be noted that Ldisc becomes larger than 
the Eddington limit in the outburst phase. Super-Eddington 
luminosities may indeed be produced during the non-steady 
evolution and/or in a non-spherical flow. In this particular 
case, this happens because of the significant transient in- 
crease in the disc temperature. The flux together with the 
disc thickness will be used in the next section to calculate 
the emergent spectra. 
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Figure 5. Same as in Fig. ^ for the effective optical depth. 
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Figure 6. Same as in Fig. ^ for the emitted flux. 

3 RADIATIVE TRANSFER IN THE DISC 
ATMOSPHERE 

In this section we outline the numerical method used in com- 
puting the evolution of the radiation spectrum emitted by 
the transonic accretion disc. As discussed previously (see 
Section ^]), we consider four stages that correspond to the 
maximal evacuation stage (model 2) and three replenish- 
ment stages (model 1, 3 and 4). 

The radial grid of each numerical disc model is made up 
of ~250 points dividing the integration domain, r^ n = 2.5 < 
r < r out — 2.8 x 10 , into a succession of comoving annular 
zones, each one 11% more massive than the one interior to 
it. The mass of the innermost zone was determined in accor- 
dance with numerical convenience. In order to compute the 
spectrum of these models, we introduce coarser radial zones 
(or rings) whose width is fixed by the requirement that their 
contribution to the total luminosity is ~ 5%. We have found 



more convenient to keep the ring boundaries coincident with 
the mesh points of the hydrodynamic calculation, avoiding 
unnecessary interpolations. The luminosity L r and average 
height h r of each ring are then computed as 



L r = ^ 27rr- i _ 1/2 Ar 4 r|.Fi„ 1/ 2 

i 

h r = ^ hj/N 



(1) 
(2) 



where n and Art are the radial mesh and grid spacing of 
the time-dependent hydrodynamic code, F i _ 1 / 2 and hi-1/2 
are the radiative flux and disc height (evaluated at the mid- 
point) and N is the number of radial points within the ring. 
Clearly, since the ring luminosity is calculated as the sum of 
the contributions from the original zones, it differs slightly 
from one annulus to another. In terms of the width Ar = 
rjv — ri and the average radius r m = (t*i + rjv)/2 of each 
ring, the luminosity can be approximately written as 



L r = 1-kt, 



lArr^^Fi/N. 



(3) 



We found that the aforementioned zoning procedure 
guarantees a sufficient accuracy when summing up all the 
rings and comparing the total disc luminosity to that com- 
puted from the output of the transfer calculation. However, 
especially in model 2, the temperature may change signif- 
icantly within a ring and this contrasts the assumption, 
which will be introduced later on, that it may be approx- 
imated as isothermal. Therefore, we further ask that the 
fractional temperature variation in each ring does not ex- 
ceed 2.5%, independently of the zone luminosity. Because 
the outermost zones of the original grid do not contribute 
appreciably to the total luminosity, they are grouped to- 
gether into one single broad ring. Typically, this procedure 
results in a total number of rings N r in 3 s = 30 — 40 (see table 
[l] for details) . 

The calculation of the radiation spectrum from each 
annulus pr oceed s i n a w ay similar to that of Shimura & 
Takahara ((1993); (1995)) for Shakura-Sunyaev discs. Each 
ring is characterized by its luminosity and height above the 
mid-plane, and is approximated as (radially) isothermal. No 
energy dissipation occurs in the disc atmosphere that is then 
in radiative energy equilibrium. All the radiative flux is pro- 
duced in the disc and enters the atmosphere from below. 
Consequently, the radiative luminosity is constant through- 
out the atmosphere, L = L r . 

The vertical structure is computed solving the cou- 
pled hydrostatic, energy and transfer equations in a plane- 
parallel, completely ionized, hydrogen layer. Solving the 
transfer of radiation in a geometrically thin plane parallel 
slab is formally very similar to solving the corresponding 
problem in a geometrically thin radial atmosphere. Thus we 
decided to model the disc atmosph ere us ing the same equa- 
tions considered by Zampieri et al. (1995) for a static, spher- 
ically symmetric, hydrogen atmosphere around a neutron 
star. The equations have been slightly modified neglecting 
general relativistic corrections, taking into account for the 
radial displacement of the ring from the center of gravity and 
suppressing the energy release within the atmosphere. The 
scattering depth r — — J K es pdz (here K es is the scat- 
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Table 1. Model parameters. 



Model 


Ldisc/LE 


Radial zones 


fave 


1 


4.8 x icr 2 


33 


1.64 


2 


1.54 


11 


1.65 


3 


1.6 x icr 2 


29 


1.65 


■1 


2.4 x icr 2 


32 


1.65 



tering opacity and p the density) is used as the indepen- 
dent variable and the inner and outer boundaries are placed 
at r in = 8 and r out — 10~ 4 , respectively. The variation of 
the height z(t) is obtained by differentiation. The radiation 
spectrum is computed from the first two moment equations 
for the monochromatic radiation energy density U v and flux 
F v , assuming that at the base of the atmosphere, where the 
T e ff > 1, the radiation field is Planckian. This set of cou- 
pled differential equations provides the run of pressure P, 
temperature T, height z, radiation energy density U v and 
flux F v as functions of depth. The density is obtained from 
the perfect gas equation of state. A thorough discussion of 
the numerical method adopted to solve these equations can 
be found in Nobili, Turolla & Zampieri (1993). 



We note that the luminosity computed by the transfer 



code is L 



r,sph 



^' Kr 'iph r G J F v dv, where r sp h is an "effec- 



tive spherical radius". Assuming J F v dv = y], Fj/N and 
L r ,s P h = L r (equation ||), we obtain 



sph = yj r m Ar/2 . 



(4) 



The spectrum emitted by the whole disc is evaluated by 
summing up the contributions from the various rings 



E 

rings 



4nr sph r G Fv . 



(•») 



4 SPECTRA FROM TRANSONIC 
ACCRETION DISCS 

Local spectra from a number of selected radial zones are 
shown in Figs. (?) and ^| for the four evolutionary stages con- 
sidered, together with the corresponding temperature pro- 
files in the disc atmosphere. We found that the emergent 
spectra are harder than the blackbody at the ring effec- 
tive temperature, T e //,r = (L r /2irr m rQAra) 1 ^ 4 (here a is 
the Stefan-Boltzmann constant). The presence of an overall 
spectral hardening is a common feature of radiative trans- 
fer in plane-parallel atmospheres and has been discussed by 
various a uthor s in different contexts (e.g. London, Howard 
& Taam ( 1986| ) for at mosp heres around X-ray bursting neu- 
tron stars; Romani (1987) and Zampieri et al. (1995) for 
atmospheres around c ooling and accreting neutron stars; 
Shimura & Takahara (1995) for the reprocessing of radia- 
tion by the top layers of Shakura-Sunyaev discs). 

The hardening is essentially caused by two distinct ef- 
fects. Models 1, 3 and 4 have an average M ~ 0.5 (in Ed- 
dington units) and are rather cool (atmospheric temperature 
Si 5x 10 6 K). Their deeper layers (r ~ Ti„) are comparatively 
hotter, T(Tjn) il T e ff. High energy photons produced there 
escape more easily, because the free-free opacity (which is 
the main source of absorption) declines rapidly with increas- 
ing frequency. The emergent spectrum is then, roughly, the 



superposition of blackbody spectra with decreasing temper- 
atures, which is peaked at m kTfcn). This mechanism is 
responsible for the spectral broadening with respect to a 
Planckian and the sensible deviation from the best-fitting 
blackbody in the local spectra of models 1, 3 and 4 (see Fig. 
|§] and the upper left panel in Fig. Q). 

The case is different for model 2, which is more luminous 
and comparatively hotter, T(Ti„) Ss 3 x 10 7 K, with an aver- 
age accretion rate M ~ 25 (in Eddington units). Local spec- 
tra at small radii are more affected by scattering and they 
appear nearly Planckian close to the maximum. The scat- 
tering depth of the inner layers (rj n ~ 10) is large enough to 
build up a Compton parameter y = 4(fcT/mc 2 ) max(r, r 2 ) ~ 
0.1 — 1. Therefore, high energy photons tend to fill the Wien 
peak at the local atmospheric temperature and photons with 
energy Jj 3fcT are sy stematically downscattered. As Shimura 
& Takahara (1995) pointed out, this explains the nearly 
Planckian shape of the spectrum close to the maximum and 
also the low energy tail in excess to the best-fitting black- 
body clearly visible in the first three spectra in Fig. [?] for 
model 2. The Compton parameter of models 1, 3 and 4 is a 
factor ~ 10 smaller, y J~ 0.1, not large enough to produce 
any spectral distortion. 

Despite the fact that local spectra may deviate appre- 
ciably from a blackbody, it is quite useful to introduce, as is 
commonly done, a hardening factor to quantify the spectral 
changes caused by radiative transfer in the disc atmosphere. 
For each radial zone the fit of the computed spectrum with a 
Planckian, cx !/ 3 /(exp hv/kT c — 1), gives the color tempera- 
ture T c . The hardening factor is defined as / = T c /T e ff and 
measures the overall energy shift of the spectrum with re- 
spect to the blackbody at T e ff, regardless of the actual spec- 
tral distortions. The run of the hardening factor in the disc is 
shown in Fig. []. Quite remarkably, we find that the harden- 
ing factor is very similar in all four evolutionary stages con- 
sidered here, and is very nearly constant through the disc. 
The typical value is / ~ 1.65 and it shows only a marginal 
i ncrea se in the outer zones. Following Shimura & Takahara 
fll995| ), we define a radially-averaged hardening factor as 



fave — 



E 

rings 



fr/Nr 



(6) 



The values of f ave for our models are listed in table |l| and the 
spectra from the whole disc are shown in Fig. [H^, together 
with the multicolor blackbody spectra at /T e // with both 
/ = fave and / = 1. As shown in Fig. [j^, with the possible 
exception of the intermediate energy range of model 2, the 
spectrum from the whole disc can be qualitatively described 
in terms of a superposition of Planckian spectra with an 
average temperature calculated using f aVB . We note that 



a similar result was found by Shimura & Takahara (199 
for Shakura-Sunyaev discs with intermediate values of the 
accretion rate and a = 0.1. 

The behavior of the hydrodynamic variables across the 
disc atmosphere is briefly illustrated here. Because the lumi- 
nosity of each radial zone is well below the Eddington limit, 
pressure scales linearly with depth. The temperature varia- 
tion across the atmosphere is rather modest: T decreases 
monotonically outwards, with the exception of the inner 
zones of model 2, where the increase is caused by the effec- 
tiveness of Compton heating/cooling in the energy balance 
(see the bottom panels in Figs. 0-^)- The atmospheric scale 
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ing the method by Cunningham (1975). We consider here 



1.55 
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Figure 9. The radial variation of the hardening factor for models 
1-4 (full, dash, dash-dotted, dash-double dotted line respectively). 



height is always quite small and the fractional increase of 
z{r out ) over the average ring height h r is only a few percent. 

Although present spectra can not be taken as represen- 
tative of the whole X-ray emission from BHCs, in that they 
do not exhibit any power-law hard tail (see section ^) , actual 
spectra as seen by an observer at infinity have been calcu- 
lated. General relativistic effects have been included, follow- 



only the case of a Schwarzschild black hole, to which the 
pseudo-Newtonian potential is the closest non-relativistic 
analog. The same calculation could be easily performed for 
a maximally rotating Kerr black hole but this seems not ap- 
propriate in the present framework. The original assumption 
of Keplerian rotation has been maintained. Deviations from 
a pure Keplerian law in our disc models are not too large 
and should not introduce major differences. The specific in- 
tensity at the top of the atmosphere (r out = 10 -4 ), needed 
to evaluate Cunningham's expression for the luminosity at 
infinity, has been obtained in terms of the first two moments 
of the radiation field as 



4>irl v (p) = cU v + 3/LtF„ = (2 + 3(j,)F v 



(7) 



where v is the frequency measured by a local observer co- 
moving with the disc and fi — cosO. Results are shown in 
Fig. hi] for different viewing angles, 8. Apart the obvious 
decrease in amplitude produced by inclination and gravita- 
tional effects, two major features are apparent in the spectra 
at infinity, particularly in the high energy tail which is pro- 
duced deeper into the potential well. First, the exponential 
cut-off is shifted to lower energies with respect to the Newto- 
nian model because of the gravitational redshift and, second, 
the spectrum tends to become harder for increasing inclina- 
tion angle. This is because the orbital motion acquires a 
larger component along the line-of-sight and this produces a 
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larger kinematical blue-shift. Finally, we note that owing to 
relativistic ray-bending, radiation from the inner zone can 
be seen also when the disc is edge-on (cos 6 = 0). 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the emission properties of 
time- dependent, slim accretion discs, calculating for the first 
time the evolution of the disc spectra along the limit insta- 
bility cycle. This work is part of an ongoing numerical pro- 
gramme to investigate the time-dependent behaviour trig- 
gered by the onset of various types of instability (thermal, 
viscou s and acoustic), as outlined in Szuszkiewicz & Miller 
(1997), and to get detailed theoretical predictions for obser- 



vational quantities relevant to accreting black holes in galac- 
tic X-ray binaries a nd AG Ns, as outlined in Szuszkiewicz, 
Turolla & Zampieri ( |200C| ). 

Away from the outburst (maximal evacuation) phase, 
the spectrum emitted by each annulus of the disc deviates 
significantly from a pure Planckian distribution, showing a 
characteristic broadening. This effect is typical of the spectra 
emitted by free- free dominated atmospheres and is caused by 
the superposition of photons emitted at different depths: be- 
cause of the energy dependence of the free-free opacity, more 
energetic photons are emitted in the deeper atmospheric 
layers, where the temperature is higher. Outburst spectra 
are more Planckian in shape around the maximum because 



the higher atmospheric temperature makes comptonization 
comparatively more effective: photons tend to form an ex- 
cess at low energies and to fill the Wien peak at high en- 
ergies. A usual way to estimate the overall hardening of a 
thermal spectrum is to compare the best-fitting Planckian 
distribution with a blackbody at the effective temperature. 
We found that, in all stages of the instability cycle, the tem- 
perature of the best-fitting blackbody is typically larger than 
the local effective temperature by a factor / ~ 1.65, indi- 
cating that the hardening factor is substantially constant 
both in radius and time. There is only a slight increase of 
/ with radius in the outer region. Therefore, although the 
mechanism responsible for the hardening is different (ther- 
mal comptonization in one case and superposition of black- 
body spectra in the others), outburst spectra have an av- 
erage hardening factor comparable to that of the spectra 
emitted during the evacuation and refilling stages. We note 
that the average value of the hardening factor for the local 
spectra in the initial state and the tendency to increase out- 
wards are consistent with what f ound for a Shakura-Sunyaev 
disc by Shimura & Takahara ((1995); / ~ 1.7) for similar 
values of the black hole mass, viscosity parameter and ac- 
cretion rate. The values of the hardening factor during the 
other stages cannot be directly compared because the struc- 
ture of the disc is rather different. 

In this investigation we focused on the properties of ra- 
diation emitted by thermally-unstable accretion discs. We 
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Figure 10. Spectra emitted by the whole disc (full lines) together with the multicolour blackbody spectra 
at T e ff (dashed lines) and at JaveT^jj (dash-dotted lines). 



remind that in the present model no energy is released in 
the disc atmosphere and that no high-energy, thermal or 
non-thermal, electron component has been included. There- 
fore, computed spectra are thermal and do not show the 
hard tail often observed in BHCs and AGNs. This must be 
taken into account in comparing present models with the 
observed spectra, in that only the evolution of the soft com- 
ponent can be reproduced. As shown in Figs, [l^ and jll], our 
model undergoes rather dramatic spectral changes along the 
instability cycle. The spectrum emitted before the instabil- 
ity sets in is typical of a standard accretion disc at moderate 
accretion rates and has a broad bump in the 0.1-1 keV band. 
During the outburst phase, the inner region becomes hotter 
and contributes most of the high energy emission in the 1— 
10 keV range: a very pronounced peak, centered around ~ 
5 keV appears in the spectrum. The spectrum then softens 
considerably during the refilling phase, up to the point that 
almost no photons are emitted above ~ 1 keV. A usual way 
to describe the observed spectra, frequently used in X-ray 
astronomy at low energy resolution, is to compute the hard- 
ness ratios. Because present spectra lack the power-law, high 
energy tail, we do not attempt to make a quantitative com- 
parison with the observed hardness ratios but simply track 
their evolution during the cycle. In Table ^ we report the 
hardness ratio HR1 (HR2) , defined as the ratio of the count 
rate in the 5-13 (13-60) keV energy band over the count 
rate in the 2-5 keV band, for the computed spectra of mod- 



els 1 and 2 (the values for model 3 and 4 are zero). The 
calculation has been performed using the RXTE response 
function and accounting for interstellar absorption (with a 
column density Nh — 2 x 10 22 cm -2 ). As shown in Table ^[ 
the hardness ratios at the beginning of the cycle are low. 
However, during outburst they increase reaching interesting 
values (HR1 ~ 1, HR2 ~ 0.03). Table | clearly indicates 
that the evolution of our model along the instability cycle is 
qualitatively very similar to the evolution o f the soft com- 
ponent in GRS1915+105 (see Belloni et al. fl200C| )). In this 
respect, it is tempting to interpret the observed soft state 
B (characterized by a higher temperature) as correspond- 
ing to the outburst spectral state of our models. The higher 
temperature is associated to the innermost part of the disc, 
accreting at much higher rates than the outer unperturbed 
part. Further investigation of this promising result is defi- 
nitely required and is presently under way. 

A few comments on the assumptions involved in our 
treatment of the accretion flow and radiative transfer should 
be made. We decided to use the slim disc approximation even 
if we are aware that realistic discs are truly 2D structures 
in which all thermodynamic quantities vary with height 
and radius. Modeling the 2D structure of accretion discs is 
much more complicated than using a ID vertically averaged 
approxim ation. However, as recently discussed by Hure & 
Galliano ( 2001 ) , the accuracy of vertically averaged models 
should be sufficient in practice for many astrophysical appli- 
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Figure 11. Observed spectra at radial infinity for different viewing angles. The full lines refer to cos = 1 
(disc face-on), whereas the dashed lines show spectra observed at cos0 = 0.75,0.5,0.25,0 (from top to 
bottom). 



Table 2. Hardness ratios for the computed spectra of models 
1 and 2. 



model 1 



model 2 



cos 9 


HRlfJ 


HR2[| 


HR1" 


HR2 6 


1 


0.021 





0.63 


0.003 


0.75 


0.042 





0.83 


0.014 


0.5 


0.055 





1.02 


0.031 


0.25 


0.064 





1.18 


0.052 





0.05 





1.31 


0.054 



HR1 = (5-10 keV)/(2-5 keV). 
HR2 = (13-60 keV)/(2-5 keV). 



cations. The slim disc model (Abramowicz et al. (198c)) re- 
tains the assumption that the disc is geometrically thin but 
the centrifugal force no longer balances gravity exactly. The 
thickness of the disc depends on both pressure and viscos- 
ity: the higher is the pressure and (in general) the viscosity, 
the thicker is the disc. In going from extremely thin discs to 
thicker ones, the effects of advection become progressively 
more important. The slim disc approach includes the effects 
of advection in a correct way while continuing to make use 
of vertical averaging. Advection is particularly important 
for the inner parts of an accretion disc around a black hole: 
while a standard Keplerian disc ends at the marginally sta- 



ble orbit, a realistic accretion flow may continue inwards and 
must cross a sonic point before reaching the event horizon. 
Finally, we note that gravity has been represented by th e 
pseudo-Newtonian potential of Paczyriski & Wiita (1980), 



which describes in a satisfactory way the dynamical effects 
of the Schwarzschild gravitational potential. The main ap- 
proximation that has been introduced in treating the vertical 
transfer of radiation in the disc concerns stationarity. This 
approximation turns out to be very good almost everywhere 
in the disc because the time-scale over which the radiation 
field changes is much shorter than the fastest dynamical time 
r/vj,, trad = (c s /cT ef f)(r/vj,) < r/v^ < r/v r , where c s , v# 
and v r are the sound, azimuthal and radial velocities. Only 
in the transonic region of the disc this approximation might 
break down. However, the radial extent of this region is usu- 
ally fairly small and does not contribute significantly to the 
total disc luminosity (^ 5% for the outburst state). The 
disc model used here neglects all the effects induced on the 
surface layers by the X-ray photons produced (or scattered) 
by an external medium. X-ray illumination by an external 
source heats the disc photosphere and leads to the formation 
of a highly ionized layer which can become Thomson thick 
for large eno ugh l uminosities (see e.g. Nayakshin, Kazanas 
& Kallman ( J2000| ); Nayakshin & Kallman j200l| )). While 
this will definitely influence the emitted spectrum shifting 
the emission at higher energies, its effect onto the hardening 
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factor are more difficult to estimate and would require a con- 
sistent calculation. We only note that present results show 
that / is not much sensitive to the disc surface temperature 
(see Table [j]). 
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